Temperature-Dependent Anomalies in the Structure of the (001) Surface of LiCu202 
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Surface corrugation functions, derived from elastic helium atom scattering (HAS) diffraction pat- 
terns at different temperatures, reveal that the Cu 2+ rows in the (001) surface of LiCu202 undergo 
an outward displacement of about 0.15 A as the surface was cooled down to 140 K. This is probably 
the first time that isolated one-dimensional magnetic ion arrays were realized, which qualifies the 
Li 1+ Cu 2+ 2 " surface as a candidate to study one-dimensional magnetism. The rising Cu + rows 
induce a surface incommensurate structural transition along the a-direction. Surface equilibrium 
analysis showed that the surface Cu 2+ ions at bulk-like positions experience a net outward force 
along the surface normal which is relieved by the displacement. Temperature-dependent changes of 
the surface phonon dispersions obtained with the aid of inelastic HAS measurements combined with 
surface lattice dynamical calculations are also reported. 

PACS numbers: 68.35.B-, 68.35.Ja, 68.35.Rh, 68.49.Bc 



Over the past decade LiCu202 has attracted consider- 
able attention because it incorporates double-chain lad- 
ders of Cu 2+ 0, which makes it a prototypical quasi- 
one-dimensional quantum spin- 1/2 magnetic system with 
competing magnetic interactions [1-4]. It still remains an 
exciting system since the recent discovery of ferroelectric- 
ity induced by an ordered helimagnetic phase, making 
LiCu202 the second cuprate to join the list of multifer- 
roics [5-8]; LiCuV0 4 [9] being the first. 

It was reported that LiCu202 exhibits a spin-singlet 
(or a gapped spin-liquid) ground state for temperatures 
T > 30 K [1]. Two successive magnetic phase transitions 
bring the system into an incommensurate long-range or- 
dered state with a helicoidal spin structure in which the 
arms of the double-chain are antiferromagnetically or- 
dered [3, 7, 8, 10-12]. The picture that emerges reveals 
the presence of a collinear sinusoidal spin-ordered phase, 
with spin polarization along the c-axis, below the tran- 
sition temperature T/v — 24.6 K. An initial modulation 
wave vector Q = (0, 0.172, 0) r.l.u. was reported, and was 
found to increase with decreasing temperature. At the 
second transition temperature Tfe — 23.0 K, the spin 
polarization acquires small components along the a and 
b axes, giving rise to the helicoidal spin-ordered structure. 
This structure is characterized by an ellipse with the heli- 
cal axis tilted by 45° from the b- axis within the ab plane, 
and with the incommensurate modulation still along the 
6-axis [10, 12]. Q continues to increase and seems to sat- 
urate at T ~ 12 K with Q = (0,0.174,0). The onset 
of the helicoidal phase induces ferroelectricity with po- 
larization along the c-axis. Moreover, recent studies by 
resonant soft x-ray magnetic scattering [8] have demon- 
strated that the long-range ordered magnetic states are 
actually 2D-like ground states, where the 2D character 
was attributed to a small but effective spin coupling along 
the c-axis, which, in turn, suppresses quantum fluctua- 
tions. These measurements also revealed the presence of 
short-range spin correlations above T/v, which eventually 



disappear at about 30 K leaving a spin liquid ground 
state at higher temperatures [8]. 

The onset of ferroelectricity below Tfe m 
LiCu202 was first reported by Park et al. [5]. Several 
models were proposed to explain the origin of electrical 
polarization in helicoidal magnetic systems. Two micro- 
scopic models were introduced to explain the mechanism 
responsible for multiferroicity in such systems, the first 
is based on the relation between the magnetoelectric 
effect and microscopic spin currents [13] and the second 
invoked the idea of inverse Dzyaloshinskii-Moriya [14]. 
In the former, the polarization is given by 

!P OC &ij x (S^ X Sj), 

where Si and Sj are the local spins at site i and site j 
respectively, and e^- is the unit vector connecting the i 
and j sites [13]. Subsequently, a phenomenological model 
based on electric polarization related symmetry- invariant 
Landau- Ginzburg terms of the form 
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was proposed [15]. In the case of LiCu202 all these mod- 
els consistently support a spiral spin with components 
in the 6c-plane, which agrees with the findings of Ref. 
[7, 8, 10-12]. 

Despite the extensive studies of bulk structural and 
magnetic properties of LiCu202 crystals cited above, to 
our knowledge, the only investigation of its surface prop- 
erties was reported by the current authors [16], and in- 
volved helium atom scattering (HAS) studies of the struc- 
ture and dynamics of the (001) surface of LiCu202 at 
room temperature. In that work, it was shown that 
the surface termination is exclusively a Li 1+ Cu 2+ 02~ 
plane. In this letter we report a surprising discovery of 
temperature-dependent anomalies observed on this sur- 
face. For temperatures below 200 K, the Cu 2+ rows were 
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found to rise above the surface plane; this, in turn, in- 
duces an incommensurate surface structure. As far as 
we know, this is probably the first time that an isolated 
one-dimensional row of magnetic ions was realized. 

LiCu202 has a layered charge-ordered orthorhombic 
crystal structure belonging to the Pnma space group. 
The primitive cell has lattice constants a=5.73, 6=2.86, 
and c=12.47 A. Single crystals of Li x Cu202 with Li con- 
tent of x ~ 0.99 ± 0.03, which can be considered as sto- 
ichiometric LiCu202 , were grown by the floating-zone 
process [17]. The experimental setup and procedures of 
the HAS measurements are described in details in Ref. 
[16, 18]. 

In this work, we present surface diffraction patterns of 
the Li 1+ Cu 2+ 02~ termination, obtained by the elastic 
HAS at temperatures below 210 K. Analysis of diffrac- 
tion intensities yielded surface corrugations, which can be 
considered as direct pictures of the geometrical arrange- 
ment of the surface atoms [19]. As we previously reported 
[16], diffraction patterns recorded between room temper- 
ature and roughly 230 K showed none of the anomalies 
discussed below. 
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FIG. 1. (a) Diffraction patterns along the a-direction ob- 
tained at T = 200, 180, 160 and 140 K. Intensities are normal- 
ized to a specular intensity of 1. Arrows indicate the presence 
of satellite peak, denoted by S, at the specified wave vectors 
q. (b) Gaussian fit (blue) to experimental peaks obtained 
at 200 K. (c) Comparison between the measured diffraction 
peaks (red) and the calculated ones (blue) . See text for details 
about the calculation. 

LiCu202 samples were cleaved in situ, under ultra- 
high- vacuum conditions at T = 200 K. Subsequently, the 
samples were cooled down in steps to T = 140 K. Fig. 
1(a) shows four diffraction patterns along the a-direction 
((lO)-direction) recorded at 200, 180, 160 and 140 K, re- 
spectively. It is clear that the intensity of the (10) first- 
order peak increases with decreasing temperature. More- 
over, a shoulder is observed on the right side of the (00) 
specular peak at 200 K and progressively separates from 
it with decreasing temperature. A Gaussian fit to the 
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FIG. 2. Surface corrugations at different temperatures with 
T decreasing. 

200 K diffraction peaks (DPs), Fig. 1(b), shows that the 
shoulder is actually a satellite peak at q = (0.11, 0) A 
Similar fits to the DPs at lower temperatures show that 
q increases monotonically with decreasing temperature 
to q = (0.25,0) A" at 140 K. Henceforth, we denote 
this satellite peak by S. Further cooling to temperatures 
below 140 K showed no change in the diffraction pattern. 

We used the hard corrugated wall model within the 
eikonal approximation to calculate the elastic scatter- 
ing intensities and followed an iterative computational 
method to fit the calculated intensities to experimental 
measurements. Details of this iterative computational 
method can be found in Ref. [16] and references therein. 
Fig. 1(c) shows a very good match between the experi- 
mental diffraction pattern at 200 K (red curve) and the 
calculated intensities of the DPs (blue rods). An impor- 
tant product of this best fit is the topology of the prim- 
itive cell manifest in the two-dimensional surface corru- 
gation function (SCF). The SCFs at T = 200 , 180, 160 
and 140 K are shown in Fig. 2, where the surface ion po- 
sitions are indicated by arrows. It is clear from the SCFs 
that the Cu 2+ ions rise monotonically with decreasing 
temperature. By calculating the difference between the 
maxima and minima of the SCFs, the largest displace- 
ment of the Cu 2+ ions is estimated to be 0.14-0.15 A 
above their bulk-like surface positions. 

The presence of the satellite peak, S, is a manifesta- 
tion of a surface incommensurate structure (SIS) along 
the a-direction. There can be two possible causes for the 
onset of the SIS: First, the appearance of surface elec- 
tric dipoles associated with the elevation of the Cu 2+ 
rows renders the surface unstable. Second, the tendency 
to maintain Cu 2+ — O 2- bond length close to its original 
value induces a lateral displacement of neighboring rows 
of O 2- ions toward the Cu 2+ rows along the a-direction. 

Fig. 3 shows five diffraction patterns along the a- 
direction sampled as the temperature was raised after 
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FIG. 3. Evolution of diffraction patterns along the a-direction 
with temperature increase above 140 K. Intensities are nor- 
malized to the specular intensity. Arrows mark satellite peak 
positions with their q values. 



cooling to 140 K. A sequence of satellite peaks emerges, 
signaling the development of the SIS phases. The loca- 
tions of these satellites are marked by arrows. The cor- 
responding wave vectors, determined by Gaussian fits, 
are recorded next to each arrow. It should be noticed 
that the satellite intensities increase with increase in 
temperature, but their wave vectors remain constant. 
The emerging SIS is characterized by satellite wave vec- 
tors: q = (0.12 ± 0.01, 0), q = (0.25 ± 0.01, 0) and 
q = (0.36 ± 0.02, 0) A . An incommensurability close 
to 9/8 was obtained by comparison with the commensu- 
rate first order wave vector of q = (1.10, 0) A , in Fig. 
1(b). All features associated with the incommensurate 
phase and the rise of the Cu 2+ rows in the diffraction 
pattern completely disappear above 220 K. 

Employing the approximate periodicity of 9/8, we car- 
ried out a peak intensity fitting procedure similar to that 
described above. The resulting SCF is shown in Fig. 
4. Moreover, we extracted the full-width-half-maximum 
(FWHM) of the specular peak from the diffraction pat- 
terns of Fig. 3; they are plotted in Fig. 5, together with 
the corresponding average domain size, as a function of 
temperature. The latter was obtained by deconvoluting 
the specular peak with an instrument transfer width [20] 
of Aki ns t = 0.07 A -1 . It is clear from Fig. 3 and 5 that 
the incommensurate domains grow with rise in tempera- 
ture. 

In Fig. 6 we plot the intensity of the ^-satellite peak at 
different temperatures, recorded during the cooling and 
heating segments of the experimental measurements. It 
is clear from Fig. 6 that the temperature dependence 
of the £-peak intensity displays a hysteresis loop with 
respect to the two segments of cooling and heating. The 
overall behavior points to a first-order commensurate- 




FIG. 4. Surface corrugation with a cross cut line for the 
incommensurate structure. 
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FIG. 5. Full-width-half-maximum (Ak) of the specular peak 
(blue), as well as the deconvoluted average domain size (red) 
as a function of temperature increase. 
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FIG. 6. The intensity of the S satellite peak. All intensities 
are normalized by setting the specular intensity to 1. The 
blue dots indicate the evolution of the corresponding g- vector 
with decreasing temperature. The red dots mark the evo- 
lution of the iS-peak intensity at the saturated q- vector, as 
obtained upon heating. The black line is a guide to eye and 
green arrows show the temperature increasing or decreasing 
procedure. 



incommensurate phase transition. 

In order to understand the experimental observations 



outlined above, we investigated the origin of the anoma- 
lous outward displacement of the Cu 2+ rows with the aid 
of detailed surface equilibrium analysis [21]. This study 
showed that at bulk-like positions the surface Cu 2+ ions 
experience a net outward force; and that this force can 
be relieved by displacing the Cu 2+ rows outward from 
their bulk positions, along the surface normal [16]. 
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tice the following trends with decrease in temperature: 
the lower branch along the 6-direction displays a gradual 
increase in frequency close to the surface Brillouin zone 
(SBZ) boundary, while the slope of the upper branch 
and its SBZ frequency gradually decrease along the a- 
direction. 

In summary, we presented experimental evidence that 
the Cu 2+ rows in the (001) surface of LiCu202 undergo 
an outward displacement of about 0.15 A as the surface is 
cooled down to 140 K; and that this displacement induces 
an incommensurate structure along the a-direction. The 
growth of the incommensurate domains was found to 
be thermally activated. The rise of the Cu 2+ rows 
was supported by surface equilibrium analysis, which 
showed that the bulk-like positions experienced outward 
forces that were relieved by the displacement. Low- lying 
surface phonon branches, associated with the motion 
of Cu 2+ and Li 1+ ions normal to the surface, exhibit 
stiffening along the 6-direction and softening along the 
a-direction as the temperature is lowered. 

M. El-Batanouny acknowledges support from the U.S. 
Department of Energy under Grant No. DE-FG02- 
85ER45222. F. C. Chou acknowledges support from Na- 
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FIG. 7. Measured inelastic events (orange dots with error 
bars) reduced to the proper surface Brillouin zone (SBZ), as a 
function of decreasing temperature, superposed on calculated 
phonon dispersion curves corresponding to different surface 
Cu 2+ displacements (Az). The calculated surface phonon 
modes are indicated by blue dots, and the gray background 
represents the projection of the bulk bands on the SBZ. Upper 
panels: a-direction. Lower panels: b-direction. 

Finally, we studied the changes in the surface phonon 
dispersion as a function of decreasing temperature. In 
these studies we employed inelastic HAS, together with 
lattice dynamical calculations based on a slab model with 
variable surface Cu 2+ positions. The procedural details 
of these studies can be found in Ref. [16, 18]. Fig. 7 
shows the results along high-symmetry directions - the 
experimental inelastic events are indicated by orange dots 
with error bars. The measurements were carried out at 
T = 200, 180, 160 and 140 K, and the surface phonon dis- 
persion curves were calculated at different surface Cu 2+ 
displacements from bulk-like positions, indicated by Az, 
which were obtained from corresponding SCFs. As was 
reported in Ref. [16], the lowest two surface phonon dis- 
persion branches (blue dots in Fig. 7) involve the motion 
of Cu 2+ and Li 1+ ions normal to the surface. We no- 
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